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Abstract 

Consider the consolidation process typical of soils, this phenomenon is expected not to exhibit a 
unique state of equilibrium, depending on the external loading and the constitutive parameters. 
Beyond the standard solution, also pore-fluid segregation, which is typically associated with 
fluidization of the granular material, can arise. Pore fluid segregation has been recognized as 
a phenomenon typical of the short time behavior of a saturated porous slab or a saturated 
porous sphere, during consolidation. In both circumstances Blot's three dimensional model 
provides time increasing values of the water pressure (and fluid mass density) at the center of 
the slab (or of the sphere), at early times, if the Lame constant fi of the skeleton is different 
from zero. This localized pore-fluid segregation is known in the literature as Mandel-Cryer 
effect. In this paper a non linear poromechanical model is formulated. The model is able to 
'Corresponding author: e-mail glulio.sclarra@uniromal.it Tel. +390644585230 Fax +390644585618 
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describe the occurrence of two states of equilibrium and the switching from one to the other by 
considering a kind of phase transition. Extending classical Biot's theory a more than quadratic 
strain energy potential is postulated, depending on the strain of the porous material and the 
variation of the fluid mass density (measured with respect to the skeleton reference volume). 
When the consolidating pressure is strong enough the existence of two distinct minima is 
proven. 
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1 Introduction 



Porous media confined into a fluid infinite reservoir, or suffering consolidation loading, can ex- 
hibit sohd-solid and sohd-fiuid phase transitions. These last can be observed because of different 



phenomena, as gravity driven solid-fluid separation of suspensions (see Burger et al. 2000) or 



sohd-fluid segregation of soils (see e.g. Nichols et al. 1994 Vardoulakis 2004a|b ). It will be the 
purpose of this paper to investigate this last transition. 



The Biot ( 1941 1955 1 three dimensional linear model describes short time pore-fluid segrega- 



tion, in case of consolidation, only for special geometries of the porous material, a slab (Mandel 



1953) or a sphere (Cryer 19631. The existence of a fluid richer stationary state of the porous 



medium can not be proven in the context of linear Biot theory endowed with Darcy solid-fluid 
viscous coupling. Recently, in-situ observations and experimental studies have pointed out the for- 



mation of compaction bands in rocks and soils (Mollema and Antonellini 1996). This phenomenon 



is typically connected with the occurrence of pore-fluid segregation (Holcomb and Olsson 2003 



Holcomb et al. 2007) in consolidation processes and, eventually, soil fluidization (Kolymbas 1998): 



under the effect of an applied external pressure, or because of gravity, solid internal remodeling can 
induce the formation of non-connected fluid-flUed cavities. Thus increasing the external loading 
causes fluid to remain trapped and therefore the fluid mass in the trapping chambers to increase 
with respect to that of the fluid flowing out of the solid matrix. The pore-fluid pressure is now 
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capable to induce imbalance of forces acting on the soil grains, so causing fluidization. Nichols 




Figure 1: Qualitative picture of the behavior of the fluid mass, m^, with respect to the consolidating 
pressure p. 



et al. ( 1994 1 have experimentally demostrated the existance of two phases in granular materials 



saturated by a fluid. A cylindrical perspex test vessel is flUed by a granular test sample and water 
is injected from the bottom through the granular layer. Tuning water pressure increases the flow 
velocity and therefore the upward drag action on the grains. At low flow velocity, when the drag 
force is smaller than gravity, the standard phase is observed: the fluid flows through the solid 
which remains undeformed. When the velocity is increased another phase shows up: the drag force 
balances gravity and fluidization of the grains occurs. 

The aim of this paper is to model solid-fluid phase transition describing in particular the 
showing up of a pore-fluid segregating state, when a consolidating external pressure is applied to 
the porous material. For this pourpose we use a simple one-dimensional model, generalizing the 
Biot theory, where the phase transition is achieved modifing the standard Biot internal energy 
functional. A suitable potential energy including the effects of the external pressure p will be 
considered. 

Also in this model the free energy depends on two flelds, the deformation e of the porous matrix 
and the density of the fluid m^, measured with respect to the solid reference volume. We shall 
describe the following phase transition: there exists pc > such that for < p < pc, that is low 



pressure, there exists a single stationary state with fluid mass m^-^^(p) and solid deformation £i{p); 
while for p > pc, that is high pressure, a second phase {m^^{p),e2{p)) appears. Fluid density rn^^ 
is greater than rn^^ and increases with the pressure p, see Figjl] 

The first solution corresponds to the case when the fluid is not confined inside the matrix and 
can flow freely back and forth from the inside to the outside of the solid. At equilibrium the density 
of the internal fluid equals that of the external infinite reservoir or, in other words the pore of the 
solid matrix are connected. The second solution, on the other hand, corresponds to the case when 
the mass density of the fluid is not that of the fluid in the external reservoir. This is due to the 
fact that the porous material starts to behave as a closed system rather than an open system, as 
the pore connecting ducts become thinner and thinner. 

This point can be supported by means of a thermodynamic argument: the flow of the liquid back 
and forth from the inside to the outside of the porous medium is a thermodynamic transformation at 
constant temperature, pressure, and volume, typical of open systems. The equilibrium is achieved 
when the internal Gibbs free energy G equals the external one; if the infinitesimal mass dm of 
fluid exits the matrix, the mass outside will vary of the amount —dm. Since the infinitesimal 
variation of the internal and the external Gibbs free energy are given by dGj = /ijdmj = ji^dm and 
dGe = /iedme = — /iedm respectively, where fi is the chemical potential of the fluid, we have that 
at equilibrium dGj = dGc implies /ij — /ie. Recalling that fi = dU /dm, with U the internal energy 
of the fluid, for reasonable choices of the function U , we have that the equality of the internal and 
external chemical potential reflects into the equality of the internal and external fluid density. 

It is quite natural that, supposed to limit our discussion to small values of the external pressure, 
the deformation of the solid is grossly proportional to the external pressure p, namely, e(p) ~ —Bp 
for some positive constant B (recall that for a compressed solid matrix the deformation is negative). 
In the unique low pressure phase, since as a consequence of the pressure p some liquid will exit the 
solid, we suppose that m^^{p) decreases proportionally to p. Concerning the second phase m^.^{p), 
the only constraint will be m^^{p) > m^-^{p). Indeed in this phase we guess the solid structure is 
modified and room for some liquid is made. 
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In spite of the appealing simplicity of these arguments, it turns out that the actual situation is 
complex to be studied. The first goal is to understand basic physical phenomena and to determine 
how phase transition in question is affected by the external pressure. 

The present paper is organized as follows. In §2 we describe the generalized Biot model and 
introduce the free energy functional. In §3, via an analytical minimization of the functional, we 
study the stationary points and their character. Finally in §4 we discuss our results. 



2 The model 

Kinematics. Let Bs, be the reference configurations of the solid and fluid components; and 
£ the Euclidean space of positions. To specify the current configuration of the system, that is 
the configuration at time t £ M, two families of diffeomorphisms, {xs.t : Bs ^ £ , t E M} and 
{4>f^ : Bs ^ B^, t E M} arc introduced. 

The current solid configuration is given by the solid placement map xs.t! for any Xg £ Bs, 
X — Xs,t(^s) is the position occupied, at time t in the Euclidean space f , by the solid material 
particle Xs- The map 4>ft^ on th^ other hand, identifies the fluid material particle X^ in B^ 
which, at time t, occupies the same current place x as the solid particle Xs- This description of the 
kinematics of the fluid is completely consistent with the Eulerean point of view adopted in standard 
fluid mechanics: the focus is not on the placement of the fluid particles, but on the particle, which 
at time t, occupies the current place x. As a consequence the reference configuration of the solid 
Bs, which in the following will be the the reference configuration of the system, is assumed to be 
a known subdomain of £, while the one of the fluid, B^, is unknown, to be determined by the 
map (p^^. As usual also the current configuration of the solid is unknown. Bearing in mind the 
definition of the map (j)^^ we shall call B :— xs.t {Bs) be the current configuration of the system. 
The so called fluid placement map Xf j '. B^ ^ £ can be constructed starting from xs,t £^nd (j)^^. 
Indeed once we set Xf ^ •= Xs,* ° "'^f ^ ^f t("''^f) represents the position occupied 



in the current configuration by the fluid particle X^. For more details we refer to Sciarra et al. 
(120081). 
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Strain. Let Fs.t '■= Vxs.t and := Vf/)^^ be the gradients of the maps xs,t and 4>^^, respectively 
(to clarify notations we remark that V indicates, in this context, spatial derivative independently 
of the domain of the map on which it operates). -Fs,t is typically named deformation of the solid. 
Since derivatives are taken with respect to Xs G Ss, those gradients are usually called Lagrangean 
gradients. We also define the gradient F^^ := Vxf^; the chain rule yields F^^{X^ — \/x{^{X^ = 
V{xsA^-^}{X^))) = ^s,t(0f;'(^f))V,^^i(Xf) = Fs^{Xs)^i,{Xsr\ where by definition X^ = 

Let Ja^t ■= \Pa.t\i with a = s,f, be the Jacobian of the transformation Xa.t measuring the 
ratio between current and reference volumes; we define the Green-Lagrange strain of the solid 
e := (Fg'^jFg t — /)/2, where / is the second order identity tensor. 

From now on we shall restrict our attention just to one-dimensional space of positions, having 
in mind to consider applications of the present model to the so-called consolidation problem (see 



Biot 1941 Terzaghi 1946 Cryer 19631. In this framework tensorial quantities restrict to scalars; 
in particular the Green-Lagrange strain reduces to e := {Jit ^ l)/2. 

Mass balance. Let po.a ■ ~^ -K, with a = s, f, be the solid and fluid reference densities. The 
total mass of each of the two components, 

M„ / go.c.{X^)dX^ (1) 

with a = s, f, is supposed to be constant with respect to time t, which means that the mass is 
conserved when passing from the reference to the current configuration. Let Qa.t, with a = s,f, 
the solid and fluid current densities, mass conservation reads 

/ go,a{Xa)dXa ^ / e'a,t(x)dx= / ga.t (Xa,t (-'^c)) ^a,t (^a )d-^a (2) 

JBc Jb Jbo, 

which in the local form, i.e. VX^ e S„, becomes ga,tiXa,t{Xa))Ja,t{Xa) — Qo.a{Xa). Note that 
the initial value q^ q of the current density of the a-th constituent can be assumed equal to the 
reference density go,a- Using the map 0^ ^ allows for introducing the solid Lagrangean mass density 



of the fluid constituent: 

mi,{Xs) go_f(0f_,(Xs))clet%^,(Xs) (3) 

The admissible deformations of the porous continuum are therefore completely known once the 
Green-Lagrange strain and the solid Lagrangean mass density of the fluid are determined. 

Overall potential. To study the equilibrium property of the system at constant temperature a 
suitable overall potential energy $, per unit volume, given by the sum of the Helmoltz free energy 
^ff and the potential of external forces, can be introduced. Since our goal is that of modeling soil 
consolidation, external loading will only consist of a pure pressure acting on the solid skeleton, 
which implies the overall potential to be $ = ^E* +pJs- 

As already noticed in a one-dimensional space of positions e = (J| — l)/2, if we restrict the 
discussion to the regime of small deformations, namely Jg w 1 and e ~ 0, we can expand around 
Js ~ 1 and get s ~ Js {geometrical linearization). 

For the sake of simplicity, from now on we shall denote with m the increment of fluid mass 
density m^, with respect to a suitable reference value For small deformations a reasonable 

expression for the dimensionless potential density of isotropic porous materials is given, in the 



framework of the Biot ( 1941 1955 ) theory, by the following quadratic form 



$^(m,e) = pe + + ^a(m - 5e)^ (4) 



where a > is the ratio between the Biot modulus M (see Coussy 2004), and the solid bulk 
modulus, while 6 > is the so-called Biot coefficient and measures the coupling between the solid 
and the fluid components; p is made dimensionless with respect to the solid bulk modulus. It is 
immediate to show that the only stationary state of Q is mg(p) = —bp and eg(p) = —p. Our 
hope is to describe the phase transition, driven by the pressure p, by considering additional third 
and fourth order terms in the overall potential. In this perspective we introduce a fourth order 
potential which reduces to the Biot one for small values of the parameters m and e, namely, in our 
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expression, the second order terms will be precisely the Biot ones. We set 



$(m,e) 



1 



am^ (3to^ - 8bem + 6&^e^ ) + * (m, e) 



(5) 
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where a — a(a, b) is a positive real function of the physical parameters a and b. While the function 
a can be chosen freely to tune the results with reasonable physical behaviors, the coefficients of the 
second order trinomial have been chosen so that ^ admits the local minimum (mi(p), ei{p)), with 
nil = bsi, describing the trivial phase of the system similar to the one obtained in the framework 
of the Biot theory. 

3 Phase transition 

We discuss the stationary states of the system which are identified with the minima of the two 
variables function m). We shall show that a point similar to the Biot one always exists, while, 
depending on the pressure, more precisely for a sufficiently large pressure, a second stationary state 
shows up. This section is devoted to the mathematical discussion of the phenomenon, its physical 
interpretation is postponed to the Section [4j 

3.1 Stationary points of the overall potential 

In order to find the stationary points of the overall potential we compute the first order partial 
derivatives of the potential energy $ 



= {iTi — be){ani^ — abem + a) 



(6) 



and 



-abm^{3be — 2m) — ab(m — be) + e + p 



(7) 
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Letting $,„ = we get the following three solutions: mi = be and 



m± = 



4a 
ab"^ 



(8) 



To get the corresponding values of e, the equation = must be solved with m = mi, m±. 
We first note that by inserting m = mi in = 0, we get 



P = -s--ab^e'^:fi{e) 



(9) 



Studying the function fi it is immediate to deduce that Q has a single real solution for any 
p > 0; we let ei such a unique solution and remark that the stationary point (wi, ei) of the overall 
potential does exist for any choice of the parameter of the model. 

To solve ([t]) with m = m±, we bound the discussion to the physically relevant region e < 0. 
Note that, if |e| > 2/{by^a/a) then 



771+ = -be 
2 



1T\ 1- 



4a 
ab-^e'^ 



< 



(10) 



which implies, in particular, that r7i_|_ > m_ > mi; remark that for \e\ large enough, m_|- approaches 
and TO_ approaches mi. By letting m — m± in Q, we get 



p = —e + ab[m±{e) — be] — ab^em'^{e) + -abmj.{e) =: f±{e) 



(11) 



Since f± are not defined around e = 0, namely, for |e| < 2/ (b^/aja), we can hope to describe some 



phase transition driven by p; indeed we guess that for p small enough the equation (111 will have no 



real solution so that (mi,ei) will be the sole stationary point of the overall potential. To explain 
our guess we draw, in Figure [2j fi, f± as function of e for the parameters specified in the captions; 
note that /_ approaches fi for \e\ large, this is consistent with the behavior of mi and m_. We 
let pc be the minimum of the function /+ and pj, :~ f+i^'c) ~ f-i^'c) with e'^, :— —2/(6-^/0/0), 
the largest value of e for which the additional stationary point of <E> appears. For p < Pc the 
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Figure 2: The graph of f± and /i is plotted. Note that the dashed hne corresponds to the graph of 
/_. The values of these functions at p determine the states of equilibrium of the porous medium. 
The used values of the constitutive parameters are a = 1/2, b — 1, a = 100. 



solutions {m± , e± ) are not real, hence the system has a single phase, the one essentially due to 
the Biot model. For pc < p < p'c the solution (m+,e+) with the good p-behavior (the one with 
smallest e+) should be interpreted as the second phase, while the other should be a saddle point 
of the overall potential. For p > p'(., the unique solution (r7i-|_,£+) should be interpreteted as the 
second phase, while (m_,e_) should be a saddle point of the overall potential. To prove that this 



interpretation is correct one should study the second order derivatives of $ (see ^3.2| 



To study the equation (11) we recall that the two functions f± are defined for e < e'q 
-2/{by^a/a) and note that 



p'c :=/±(£c) = 



—^(2 + lbW >0 
by/a/a ^ o 



(12) 



To study the e — >■ — oo limit, we recall that \/l — x = 1 — x/2 — + 0{x^), and using (10 1 
we get 

1. r.. , 2a , 2a- ^^^^ 



m±{e) = ^be 



where 5-,± is the Kronecker delta, such that (J-,- = 1 and — 0. Substituting those expansions 



10 



in the definition (11) of /± and keeping track of the not vanishing terms we get 



/+(e) - -e(l + ab^) and /_(e) - -e - -ab'^e^ 



(14) 



Note that for e — > — oo the function /_ behaves precisely as /i 



We now compute the first order derivative. Recalhng the definition (11) of f±, see equation 
([6]), and the fact that m± are solutions of the equation $,„ = 0, we get 



/4(e) = - 1 — ab'^ — ab'^m^ (e) + bm'^ (e) [a — 2abem± [e] + 2am\ (e)] 



(15) 



Note the the first three terms are clearly negative, so that the sign of the derivatives depends 



essentially on the sign of the fourth term. We then recall (10), which is valid for e < 0, and 
compute 



m±(£) = -m±{£) =F —r-i^ 1 



4a 



y.b'^e'^ 



-1/2 



(16) 



By using (16) and (15) it follows immediately that /±(e) -> ±oo for e e'q from the left. 



Notice that m!_{e), for e < 0, is positive. Hence, by (15) it follows that f'_{s) < for e < e'^ 



The function /_ decreases monotonously from +oo to pj, as e goes from — oo to e'^. We then 
conclude that the stationary point (m_, £_) of the overall potential exists and is unique for p > Pq. 



The study of is slightly more difficult. First of all we recall that > m+{e) > be, see (10) 



hence, using the definition ( 11 ), we get 



/+(e) = — e + a6(m+(e) — be) + abm'^{e) -m+{e) ~ be 



> 



(17) 



Thus /+ (e) is a positive function tending to +oo for e — > — oo and approaching /+ (sc) > 0; see ( 12 ) , 
with positive slope. Hence there must be at least a minimum of the function /+(e) in the region 
e E (— oOjEcJ. Moreover, it is possible easy to see that the equation f'^_{e) = is biquadratic in e. 
Hence it can have either zero or one or two negative solutions, the only possible case, compatible 



11 



with the above mentioned properties of /+, is that the negative solution is unique. 

Finally, we set pc '■= /+(£c) > and remark that for p < pc the equations /±(e) = p have no 
real solution. For pc < p < p^, the equation /-(e) = p has no real solution, while /+(e) = p has two 
real solutions {m\,e\_) and (m^,e^) with e\ < e\. For p > p'^, both /-(e) = p and /+(e) = p 
have a unique real solution, respectively denoted by (m_,e_) and {m^,e\). According to the 



above depicted scenario, see the discussion below (111, we expect that {m^,e^) is a minimum of 
the overall potential, while (TO^,e^) and (m-,e-) are saddle points. The validity of this guess 
will be proven in the next section. 

3.2 Character of the stationary points 

We study, now, the character of the stationary points of the overall potential. Computing the 
second order derivatives of $ with respect to m and e we get 

^mm = {am^ — abem + a) + {m — be)(2am — abe) 

= -ab + 2abm{be - m) (18) 
= 1 + ah^ + ab'^m? > 

Using that mi = feei, we are able to compute the Hessian H{mi,ei) = a(l + am\b^) > and, 
thus, conclude that (mi,ei) is a local minimum of Hence it represents a stationary state of 
the model. 

We have to study, now, the properties of the stationary points {m±,e±). In view of this we 



give a nice expression of the Hessian computed in (TO±(e),e). By using (181 and recalling that 
m±(e) is obtained by solving the equation = 0, see (|6]), we get 

^mm{m±{e),e) = a{m±{e) - be){2m±{e) - be) 

^me{m±{e),e) = ab (19) 
$£e(m±(e),e) = 1 + b'^{a + a{m±{e))'^) > 



12 



so that 



T~L±{^) ■= 'H{m±{e),e) 



By using we get 



[1 + 52(a + a(TO±(e))^)][a(rn±(e) - be){2m±{e) - be)] ~ a^b"^ 
[1 + b"^ abem±{£))][a{'m±{£) — 6e)(2m±(e) — foe)] — o^b"^ 



(20) 



-H±(£) = -2a - aH"^ + ab^e^ ab^ [1/2 + ab^)£^/7^^^^1ilij{^ 



(21) 



Since we have focused our discussion on the case e < 0, we have that 



U±{e) = -2a - aH'^ + ab^e^/2 ± ab^ {1/2 + ab^)e^ ^^1 - Aa/{ab^e^) 



(22) 



We use the above expression (22) to prove that the stationary point (m_,e_) is a saddle point. 



First we note that if such a stationary point exists then it must be necessarily £_ < £c, so that 
4a/(afe^e^) < 1. Remarking that, for < a; < 1, one has \/l — x > 1 — x, by (22) we get the 
bound 



•H_(e_) < -2a-a^b^ + ab^el/2-ab^{l/2 + ab'^)e'i[l-4a/{ab'^e'i)] 
= ab^[3a - ab^el] < ab^ (3a - 4a) = -a%^ < 

Where the last bound follows from the fact that, since £_ < e'q, we have ab^e^ > 4a. 



(23) 



In order to study the character of the stationary points (m*^,£*|_), with i — 1,2, we study the 
fmiction H+ie) for e < 0. Recall (22) and note that H+{e'(.) = —a^b^ < and H+ie) — > +oo for 



e — )■ —oo. By computing the first derivative of 'H+{e) with respect to e we get 



(24) 



which is clearly negative for e < e'(.. Thus, we have that 'H+(e) decreases from +oo to —a b as e 
goes from — oo to e^. Denoted by e the unique negative zero of 'H^{e), we prove that £ = ec, that 
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is £ is the minimum of the function f+{s). Note, indeed, that by definition 



H+{e) = $„„(m+(£),£)$ee("i+(e),e) " *™^(m+(£),£) (25) 
Recall that the function m^^(e) is imphcitely defined by the equation <I>m(m_(_(£), £) — and note 



that (11) is equivalent to f+{e) = — <i>e(rn+(£), £) + p. By using the chain rule and the implicit 



function theorem we get that 

H+{e) 



*mm(m+ (£),£) (^^Q) 



where in the last step we have used ( 25 1 



The above equality ensures that £ = £c. Finally, since, £^ < £c < £+, we have that {m^,e\_) 
is a minimum of ([5| while (m^,£^) is a saddle. 



4 Results 

In this section a discussion on the behavior of the minima of the overall potential ^ , when varying 
the external pressure, will be presented. A parametric analysis in dependence of the coefficients 
introduced in the Biot model, a and b, as well as of the additional coefficient a (a, b) multiplying 
the fourth order terms of ^ will be also developed. 

According to the experimental results on fluidization of soils, cited in the introduction (see 
Nichols et aTj |1994[ [Holcomb and 01sson[ |2003| [Holcomb et~aH |2007| [Vardoulakis] |2004a|b[ ), 



the proposed model is capable to exhibit, in the presence of an applied consolidating external 
pressure, an additional stationary state in the (£, m) plane beyond the one corresponding to classical 
consolidation. Increasing of the external pressure, acting on the porous medium, induces, on one 
hand, part of the fluid to flow out of the skeleton, on the other, part of the fluid to be segregated 
into not connected cavities of the solid matrix, when the external pressure overwhelms a critical 
value. 
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The arising of the additional state of equilibrium, associated with pore-fluid segregation, is 
explicitly depicted in Fig. ([s]), where the opposite of the total stress a = d"^ /de and the fluid 
chemical potential /i = d'^ /dm are plotted against the kinematical parameters e and to. The 
projection over the (e, to) plane of the curves obtained cutting the —a surface with the horizontal 
plane at p and the /i surface with the (e, m) plane itself have more than one mutual intersection, 
each of them identifies one of the stationary points of the overall potential ^ . In particular the 
dotted (dashed) lines in Fig. ([s]) (Fig. Q) correspond to the solutions of = 0, on the other 
hand the solid ones correspond to the solution of $j = 0, for different values of p. 




Figure 3: Cutting the —a surface with the horizontal plane at p identifies the root locus of 4>e = 0; 
in particular the two solid lines identify the root loci associated with the critical value of the applied 
pressure and that associated to a larger value of p. The dotted lines describe the root locus of 
$,„ = 0, obtained cutting the /i surface with the (£,m) plane. 

As already noticed, the critical pressure pc is the smallest value of the consolidating loading for 
which the pore-fluid segregation solution shows up. 

With another point of view, the stationarity conditions stated by eq. ^ and ((7| describe two 
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Figure 4: This is a detail of the projection on the (s, m) plane of the intersection of —a and /z 
with the horizzontal plane at p and zero, respectively. In particular (e^,m^), i — 1,2, identify the 
two stationary points of ([s]) when p = pc, while (ef^, mf^), i = 1,2, 3, identify the three stationary 
points of ([5]) when p > pc (super-critical conditions). 
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Figure 5: The dotted lines correspond to the intersection between —a and the root locus of — 0, 
the solid lines are the projections of these curves on the stress-strain plane, i.e. the /i and f± 
curves. 

different surfaces in the (£,m) plane, the intersection of which gives rise to a curve in the three- 
dimensional space, the dotted curve in Fig. ([5]), playing the role of the curve of the fluid pressure 
in terms of the specific volume, in the case of Van der Waals' model of liquid-vapor coexistence. 
It is worth to notice that the projection of these curves on the (— cr, e) plane provides the plot of 
the /i and f± functions, depicted in Fig. As in the case of Van dcr Waals' model sectioning 
this curve in the three dimensional space with an horizzontal plane at p identifies one stationary 
(equilibrium) point, if p < pc, or three stationary points, if p > Pc, two of which corresponds to 
classical consolidation and pore-fluid segregation solutions, respectively. 

Fig. ^ explicitly shows the parametrization with pressure of the states of equilibrium: when 
the applied external pressure reaches the critical value Pc a new state of equilibrium appears in 
which the variation of the fluid mass and the negative strain increase. This is what is generally 
called fluid-segregation: under the applied external pressure the porous material undergoes to 
consolidation, however part of the fluid remains trapped into non-connected cavities of the skeleton, 
the formation of which is possibly due to the thinnering of duct connections of the pore-space. 

Note that all the pictures have been drawn for the parameters specifled in the caption of Fig. ([2]). 
However it is important to stress that, according to the general discussion of SjSj no variation in 
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Figure 6: The graph of /i corresponds to the classical behavior of consolidating soils, the neg- 
ative strain progressively increases and the mass of the fluid progressively decreases when the 
external pressure raises up; the graph of /+, on the other hand, describes the phenomenon of 
fluid-segregation, the dashed line (graph of /_) corresponds to the saddle point of the overall 
potential. 

the main features of the considered phase transition, from the classical consolidating soil phase 
towards the pore-^fluid segregated one, can arise when tuning the constitutive parameters a, h and 
a. The only constraint on the showing up of the new equilibrium being the impenetrability of the 
skeleton (e > — 1/2). It is worth to notice that the absolute values of both the free energy of the 
pure fluid and that associated with solid-fluid coupling in the overall free energy ([s]) 

*g-(e, m) := l/4,am'^ + l/2am^ and ^I'gf (e,m) -2/3a&m^e + l/2a6^m^e^ - ahme (27) 

progressively decrease for increasing values of the applied external pressure, along the pore-fluid 
segregating equilibrium path, see Fig. ([7|. Apparently this is due to the increasing incapability 
of the pore fluid to flow out of the skeleton. Conversely, the free energy of the skeleton, 5'ss := 
1/2 (l + a6^) along the pore-fluid segregating equilibrium is larger than the one associated to 
classical consolidation, the trapped fluid inducing a virtual stiffening of the solid matrix. 

Phase transition is therefore completely characterized in terms of the dependence of the critical 
pressure on the ratio between the Biot modulus and the bulk modulus of the skeleton, a, and 
the Biot coefRcient b. In the frame of physically meaningful assumptions on these parameters 
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Figure 7: The solid lines indicate the pure solid, the pure fluid and the coupling energies associated 
with the soil consolidating equilibrium path; the dashed ones indicate the same energy contributions 
associated with the pore-fluid segregation equilibrium. 
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Figure 8: Picture of the variation of the critical pressure with the Biot coefficient b, parametrized 
by a; both a and b range in the open interval (0, 1). 
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(a < 1 and 6 < 1) the critical pressure pc exhibits the behavior depicted in Fig. (|8|. Increasing a, 
keeping h fixed, impUes the critical pressure to increase and therefore the phase transition to occur 
for higher values of the consolidating pressure. Indeed increasing a means increasing of the bulk 
modulus of the fluid; the inverse of Blot's modulus is an affine function of the inverse of the fluid 
bulk modulus, with respect to that of the solid. The fluid becomes more stiff and therefore does 
not allow the applied external pressure to force shrinkage of the ducts inside the porous medium; 
as a consequence flowing of the fluid out of the solid matrix is enhanced. 

Increasing 6, keeping a fixed, on the other hand, implies the critical pressure to descrease and 
therefore the phase transition to occur for lower levels of the consolidating loading. This is due to 
the fact that increasing h enhances solid-fluid coupling. 

Apparently no physical reasoning can be developed in order to discuss the behavior of the 
system when tuning a, however it is clear that only for sufficiently high values of this parameter 
phase transition can be appreciated. 

References 

M. Biot. General theory of three dimensional consolidation. Journal oj Applied Physics, 12 (2): 
155-164, 1941. 

M. Biot. Theory of elasticity and consolidation for a porous anisotropic solid. Journal of Applied 
Physics, 26 (2):182-185, 1955. 

R. Biirger, W.L. Wcndland, and F. Concha. A mathematical model for sedimentation-consolidation 
processes. Zeitschrift fiir Angewandte Mathematik und Mechanik, 80:79-92, 2000. 

O. Coussy. Poromechanics. John Wiley & Sons, Ltd, 2004. 

C.W. Cryer. A comparison of the three-dimensional consolidation theories of Biot and Terzaghi. 
QuaH. J. Meek, and Applied Math., 16 (4):401-412, 1963. 



20 



D. Holcomb and W.A. Olsson. Compaction localization and fluid flow. Journal of Geophysical 
Research, 108 (B6):2290, 2003. 

D. Holcomb, J. W. Rudnicki, K. A. Issen, and K. Sternlof. Compaction localization in the earth 
and the laboratory: State of the research and research directions. Acta Geotechnica, 2 (1):1-15, 
2007. 

D. Kolymbas. Behaviour of liquified sand. Philosofical Transactions of the Royal Society of London 
A, 356:2609-2622, 1998. 

J. Mandel. Consolidation des sols (etude mathematique). Geotechnique, 3:287-299, 1953. 

P. N. Mollema and M. A. AntoneUini. Compaction bands: a structural analog for anti-mode 1 
cracks in aeolian sandstone. Tectonophysics, 267:209228, 1996. 

R.J. Nichols, R.S.J. Sparks, and C.J.N. Wilson. Experimental studies of the fiuidization of layered 
sediments and the formation of fluid escape structures. Sedimentology, 41:233-253, 1994. 

G. Sciarra, F. dell'Isola, N. laniro, and A. Madeo. A variational deduction of second gradient 
poroelasticity: Part I: General theory. Journal of Mechanics of Materials and Structures, 3 (3): 
507-526, 2008. 

K. Terzaghi. Theoretical Soil Mechanics. John Wiley & Sons, Ltd, 1946. 

I. Vardoulakis. Fiuidization in artesian flow conditions: hydromechanicafly stable granular media. 
Geotechnique, 54:117-130, 2004a. 

I. Vardoulakis. Fiuidization in artesian flow conditions: hydromechanicafly unstable granular 
media. Geotechnique, 54:165-177, 2004b. 



21 



